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Abstract 

An  approximate  non-asymptotic  theory  of  weak 
normal  shock-unseparated  turbulent  boundary  layer 
interactions  is  given  which  includes  the  effect 
of  surface  mass  transfer.  The  results  of  a para- 
metric study  of  Reynolds  number  effects  on  various 
interaction  properties  without  mass  transfer,  in- 
cluding skin  friction,  are  first  presented  along 
with  detailed  comparisons  with  experiment  in- 
cluding supercritical  airfoil  data.  The  extension 
and  application  of  the  theory  to  include  moderate 
amounts  of  either  suction  or  blowing  through  the 
surface  is  then  discussed,  especially  as  regards 
the  influence  of  mass  transfer  on  the  skin  friction 
behavior.  As  a consequence  of  its  influence  on 
the  boundary  layer  profile  shape  away  from  the 
wall,  small  suction  (-mw/i’eue  < 5 x 10-4)  appreci- 
ably reduces  the  upstream  influence  and  thickening 
effects  of  the  interaction  but  hastens  the  onset 
of  separation  in  the  shock  foot  region. 


Nomenclature 

B Mass  transfer  parameter,  m w/Pe0ue 

Cf  Skin  friction  coefficient, 

V/^ef 

L Distance  to  undisturbed  shock  location 

m Mass  flux  rate  (per  sec.)  across  sur- 

face  (=  w 
0 0 

M Mach  number 

aP  Pressure  rise  across  incident  normal 

shock 

p'  Interaction  pressure  perturbation 

(P  - P,) 

Re. ,Re  Reynolds  numbers  p„u_L/p  , p„u  6/p 

l 6 e*  e0  eo  eo  eo 

T Absolute  static  temperature 

u,v  Flow  velocity  components  along  x,y 

respectively 

x,y  Streamwise  and  normal  coordinates, 

respecti vely 

y Ratio  of  specific  heats 

5 Undisturbed  boundary  layer  thickness 

6*  Displacement  thickness 

•1  y/« 

An  Interactive  displacement  growth, 

A6*/«0* 

u Coefficient  of  viscosity 

v Kinematic  viscosity  coefficient,  u/p 

p Density 

t Shear  stress 

e*  Momentum  thickness 

Subscripts 

1,2,3  Interaction  regions.  Fig.  2 

e Edge  of  boundary  layer 

0 Undisturbed  solid  wall  boundary  layer 

conditions 


* Professor  of  Aerospace  and  Ocean  Engineerin 


Subscripts  (Continued) 

ref  Eckert  reference  temperature  value 

w Wall  surface  conditions 


1.  Introduction 

The  study  of  transonic  shock-turbulent  bound- 
ary layer  interactions  is  important  in  the  aero- 
dynamic design  of  high  speed  aircraft  wings,  tur- 
bine and  cascade  blades  in  turbomachinery,  and 
airbreathing  engine  inlets  and  diffusors.  Con- 
sequently, the  control  and  suppression  of  inter- 
action effects  in  these  applications  by  suction 
or  blowing  is  of  interest  since  boundary  layer 
control  (BLC)  by  surface  mass  transfer  has  achieved 
practical  status.  Experimental  studies  have  estab- 
lished the  value  of  using  distributed  suction 
normal  to  the  surface  for  shock-boundary  layer 
interaction  separation  suppression  on  wings' *2, 
porous  wind  tunnel  walls®  and  within  supersonic 
inlets^.  There  is  also  interest  in  normal  blowing 
effects  in  connection  with  transonic  shock-bound- 
ary layer  interactions  on  wings®  and  on  mass  trans- 
fer - cooled  hot  turbine  blades®.  Moreover,  the 
study  of  interactions  in  the  presence  of  dis- 
tributed normal  mass  transfer  is  of  fundamental 
interest  in  its  own  right?  and  provides  an  ideal- 
ized model  for  optimum  performance  estimates  of 
BLC  systems  on  wings  and  flaps  and  rational  cri- 
teria for  separation  prevention  by  suction®. 

Although  some  basic  theory  for  the  transonic 
shock-turbulent  boundary  layer  interaction  prob- 
lem has  accumulated  for  the  zero  mass  transfer 
case  only  rather  crude  overall  integral 

methods  are  presently  available  to  treat  the 
effects  of  suction  or  blowing**;  a more  detailed 
theory  including  mass  transfer  is  thus  desire- 
able  to  provide  a proper  analytical  framework  and 
understanding  in  the  aforementioned  applications. 
Likewise,  there  has  evidently  been  no  systematic 
study  of  attendant  Reynolds  number  effects  on  un- 
separated interactions  at  realistic  values  per- 
taining to  practical  applications,  although  these 
effects  may  be  significant.  The  purpose  of  this 
paper  is  to  present  the  results  of  a combined 
study  of  these  two  features  for  the  case  of  non- 
separating flow  with  normal  wall  (unvectored)  mass 
transfer. 

Our  approach  is  based  on  extending  a pre- 
viously-developed approximate  non-asymptotic 
theory  of  weak  normal  shock-turbulent  boundary 
layer  interactions*1  on  the  premise  that,  not- 
withstanding the  existance  of  powerful  numerical 
methods,  there  will  be  a continuing  need  for  such 
analytical  methods  that  delineate  the  es>ential 
physical  features  and  parametric  trends  of  the 
problem  at  realistic  Reynolds  numbers.  In  Sec- 
tion 2 the  basic  formulation  and  features  of  the 
theoretical  model  Including  surface  mass  transfer 
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effects  are  given.  Section  3 then  discusses  typi- 
cal numerical  results:  first  for  zero  mass  trans- 
fer where  heretofore  - unpublished  Reynolds  num- 
ber effects  including  those  on  upstream  influence 
and  skin  friction  are  presented  plus  several 
comparisons  with  experimental  data;  results  show- 
ing the  mass  transfer  effects  on  important  inter- 
action properties  are  then  given.  Section  4 con- 
cludes with  a discussion  of  the  limitations  of 
the  theory  and  recommendations  for  further  studies. 


2.  Theoretical  Formulation 
2.1  Basic  Features  of  the  In teraction  Flow  Model 

It  is  well-known  experimentally  that  when 
separation  occurs,  the  disturbance  flow  pattern 
associated  with  normal  shock-boundary  layer 
interaction  is  a very  complicated  one  involving  a 
bifurcated  shock  pattern1?,  whereas  the  unsepara- 
ted case  pertaining  to  turbulent  boundary  layers 
up  to  M]  < 1.3  has  instead  a much  simpler  type  of 
interaction  pattern  which  is  more  amenable  to 
analytical  treatment  (see  Fig.  1).  With  some 
judicious  simplifications,  it  is  possible  to  con- 
struct a fundamentally-based  approximate  analytical 
theory  of  the  problem  in  this  latter  case.  For 
the  sake  of  orientation  and  completeness,  a brief 
summary  of  this  theory  will  now  be  given  (full  de- 
tails can  be  found  in  Ref.  11). 

The  flow  consists  of  a known  incoming  iso- 
baric  turbulent  boundary  layer  profile  M0(y)  sub- 
jected to  small  transonic  perturbations  due  to  an 
impinging  weak  normal  shock.  In  the  practical 
Reynolds  number  range  of  interest  here  [Re^ 

0(106)]  we  purposely  employ  a non-asymptotic 
disturbance  flow  model  in  the  turbulent  boundary 
layer  patterned  after  the  Lighthill-Stratford- 
Konda  double-deck  approach! 3-1 5 that  has  proven 
highly  successful  in  treating  a variety  of  other 
problems  involving  turbulent  boundary  layer  re- 
sponse to  strong  rapid  adverse  pressure  gradi- 
ents1^'1?. This  was  done  because  of  the  large 
body  of  turbulent  boundary  layer-shock  interaction 
data  that  strongly  supports  such  a model  in  this 
Reynolds  number  range!8-28  (including  the  transon- 
ic regirie?9-31)  and  because  of  the  findings  of  a 
separate  general  theoretical  study32  showing  that 
asymptotic  theory  results  for  very  high  Reynolds 
numbers  5,1 j.  although  rigorous  in  this  limit,  do 
not  extrapolate  down  to  the  present  Reynolds  num- 
ber range.  The  resulting  flow  model  consists  of 
an  inviscid  disturbance  flow  surrounding  a non- 
linear shock  discontinuity  and  underlaid  by  a thin 
viscous  disturbance  sublayer  as  schematically  il- 
lustrated in  Fig.  2.  The  introduction  of  some 
further  simplifications  (including  the  assumption 
of  small  linearized  disturbances  ahead  of  and  be- 
hind the  nonl inear  shock  jump*  plus  neglect  of  the 
detailed  shock  structure  within  the  boundary  lay- 
ers, which  give  accurate  results  for  the  overall 


*As  far  as  the  overall  interaction  solution  for 
1 05  < Rej_  <_  1 08  is  concerned,  these  nonlinear 
shock  jump  conditions  plus  the  various  non- 
uniform  viscous  flow  effects  within  the  boundary 

layer  reduce  the  lower  Mach  No.  limit  other- 
wise pertaining  to  the  linearized  supersonic 
theory  in  purely  inviscid  potential  uniform  flow- 
see  the  AIAA  paper  version  of  Ref.  11  for  more 
detailed  discussion. 


properties  of  engineering  interest  provided  Mi  is 
not  too  close  the  unity  (M-j  ; 1.05),  then  yields 
an  approximate  analytical  solution  consisting  of 
linearized  potential  supersonic  disturbance  flow 
in  region  1 plus  the  subsonic  disturbance  flow  in 
quadrant  3 caused  by  the  interaction-generated 
interface  displacement  i3(x)  and  the  postshock 
perturbations  along  x - 0+  due  to  the  impingement 
of  region  1 Mach  wave  disturbances  on  the  shock; 
beneath  this  is  a rotational  inviscid  disturbance 
flow  region  2 inside  the  boundary  layer  under- 
laid by  a thin  shear  stress-disturbance  sublayer 
which  carries  the  major  upstream  influence,  dis- 
placement thickness  effect  and  skin  friction  per- 
turbation of  the  interaction.  In  the  Reynolds 
number  range  of  interest  here,  we  can  approximate 
this  thin  sublayer  to  lie  within  the  linear  part 
of  the  basic  velocity  profile  where  U0  = (iw/>‘w )y 
(i.e.,  within  the  laminar  sublayer  of  the  tur- 
bulent profile**);  consequently,  a known  solution^ 
for  this  sublayer  can  be  extended  to  the  present 
problem  to  obtain  the  displacement  effect  (effec- 
tive wall  position  seen  by  the  overlying  inviscid 
disturbance  flow)  and  the  corresponding  disturb- 
ance skin  friction  caused  by  the  pressure  dis- 
turbance field. 

The  matching  of  these  regional  solutions 
yields  linear  integral  equations  that  can  be 
readily  solved  by  operational  methods11  for  the 
disturbance  pressure  along  both  the  boundary 
layer  edge  and  wall.  The  remaining  interactive 
flow  properties  can  then  be  determined  in  terms 
of  pw' , including  the  interactive  growth  in  dis- 
placement thickness  and  the  corresponding  skin 
friction  perturbation.  It  is  noted  that  this 
latter  perturbation  solution  has  been  recently 
extended  to  include  the  region  downstream  as  well 
as  upstream  of  the  shock  and  to  include  non- 
linear inertia  effects  in  an  adverse  pressure 
gradient  using  the  general  non-dimensional  wall 
shear-pressure  solution  ahead  of  separation  given 
by  triple  deck  theory36  (converted  to  turbulent 
flow  by  expressing  all  results  in  terms  of  Cf0  in- 
stead of  ReL);  the  results  yield  the  approximate 


(1) 


where  the  non-dimensional  function  F is  essential- 
ly unity  ahead  of  the  shock  x <0  and  then  de- 
creases behind  it,decaying  slowly  like  F ■*» 
(x/60)*1/3  far  downstream.  Eq.  1 indicates  that, 
depending  on  Reynolds  number  (Cf0),  a sufficiently 
strong  interactive  pressure  rise  can  cause  in- 
cipient separation  [Cf(x)  - 0]  near  the  shock. 

This  solution  contains  the  essential  global 
features  of  the  mixed  transonic  character  of  the 
non-separating  normal  shock-turbulent  boundary 
layer  interaction  problem  including  the  signifi- 
cant lateral  pressure  gradient  effects,  upstream 


**As  shown  below,  this  gives  good  results  for 
Ret  ~ 10f  - 10a;  however,  it  begins  to  signifi- 
cantly break  down  for  ReL  > 108  where  the  effect 
of  the  logarithmic  portion  of  the  profile  on  the 
sublayer  solution  becomes  significant32. 
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influence  and  interactive  skin  friction  for  an 
arbitrary  input  turbulent  boundary  layer  profile. 
Moreover,  many  detailed  comparisons  with  all 
available  experimental  data  (Ref.  31  and  below) 
plus  recent  streamlining  of  the  computation  pro- 
gramed", have  shown  the  method  to  give  a very 
good  account  of  all  the  important  engineying 
features  of  the  interaction  for  ReL  < 10®  at  low 
computational  cost.  Hence  the  theory  provides  a 
sound  basis  for  interpreting  experimental  data  on 
unseparated  flows  and  for  further  extensions,  in 
particular  to  allow  mass  transfer  through  the  sur- 
face. 

2.2  The  influence  of  Mass  Transfer 

The  influence  of  surface  mass  transfer  on  the 
interaction  solution  is  two-fold:  (1)  it  alters 
the  incoming  undisturbed  flow  (on  which  the  dis- 
turbance solution  depends)  by  changing  iw  and 
the  profile  shape  away  from  the  wall;  (2)  it  fur- 
ther introduces  new  mass  transfer- induced  terms  in 
the  disturbance  equations.  Now  the  secondary  ef- 
fects (2)  may  in  fact  be  neglected  with  good  ap- 
proximation under  the  assumed  conditions  of  small- 
to-moderate  normal  mass  transfer  rates  ( bi*/.  e u^|- 
1 0- 3 ) typical  of  practical  applications,  according 
to  the  following  considerations.  Under  the  con- 
tinued assumption  that  the  viscous  disturbance  sub- 
layer lies  within  the  laminar  sublayer  region  of  a 
turbulent  boundary  layer,  the  mass  transfer  effect 
in  the  leading  approximation  does  not  introduce  any 
curvature  into  the  Uo(y)  profile  but  only  alters 
its  slope  (ij;  consequently  in  the  viscous  dis- 
turbance sublayer  perturbation  equations  those 
terms  proportional  to  d2l^,/dy2  and  vWo  (which  in 
turn  is  proportional  to  d;>Uo/dy2  via  the 
mean  fir'*  mentum  equation  near  the  wall)  can  be 
neglectcu,  leaving  the  form  of  these  equations  un- 
changed. Likewise,  under  the  continued  assump- 
tion that  turbulent  fluctuations  are  uncorrelated 
with  the  interactive  disturbances,  the  explicit 
new  terms  in  the  overlying  rotational-inviscid 
perturbation  equations  that  are  proportional  to 
vWo  can  be  neglected  also,  since  detailed  studies 
of  the  hydrodynamic  stability  equations-’"  (which 
are  very  similar  to  those  of  the  present  problem) 
have  shown  that  these  terms  have  an  altogether 
negligible  effect  throughout  a high  Reynolds  num- 
ber parallel  shear  flow  boundary  layer  unless  the 
surface  mass  transfer  is  quite  large  (approaching 
blow-off).  Thus  to  a consistent  degree  of  first 
approximation,  the  form  of  both  the  viscous  and 
inviscid  disturbance  equations  is  unchanged  by 
moderate  blowing  or  suction  provided  their  effect 
on  the  mean  flow-based  coefficients  in  these  equa- 
tions included.  It  is  reemphasized  that  the 
primary  mass  transfer  effects  on  the  viscous  sub- 
layer field  and  its  thickness  are  thus  fully 
accounted  for. 

The  present  interactive  perturbation  solution 
in  principle  may  be  used  with  any  mean  turbulent 
boundary  layer  profile  input  and  hence  could  be 
coupled  with  either  an  experimental  measurement 
or  any  desired  state-of-the-art  numerical  pre- 
diction code.  Here,  to  bring  out  clearly  and 
efficiently  where  the  various  mass  transfer  ef- 
fects enter,  we  have  chosen  an  accurate  analytical 
profile  model4"  that  has  proven  especially  well- 
suited  to  such  non-uniform  flow  perturbation  prob- 
lems^'. We  assume  for  simplicity  that  the  blowing 
or  suction  is  on  the  average  uniform  and  normal  to 
the  wall,  that  its  streamwise  extent  is  large  com- 


pared to  the  short  interaction  range,  and  that  it 
extends  far  enough  upstream  to  have  established  a 
well-defined  local  equilibrium  profile  in  the  in- 
coming boundary  layer.  Then  the  mass  transfer 
effect  on  turbulent  skin  friction  in  terms  of  the 
basic  parameter  B - ii>w/ can  "e  described  by 
the  relation4? 


where  Tref/Te  - 1 + 0.038Me  + 0.50(TW/Te  - 1)  for 
> = 1.4,  and  Cf  is  the  zero-blowing  value  (we 
used  a reference  temperature-based  Schulz-Grunow 
relation’).  According  to  Eq.  (2)  suction,  for  ex- 
ample, increases  the  skin  friction. 

The  corresponding  mass  transfer  effect  on  *0 
can  be  estimated  as  follows.  Since  it  is  well- 
known  that  the  momentum  thickness  to  boundary 
layer  thickness  ratio  and  the  Crocco  energy  equa- 
tion solution  are  both  insensitive  to  moderate 
amounts  of  mass  transfer^,  We  can  use  the  follow- 
ing approximate  zero  blowing  relationship  based 
on  a power-law  (U0  ***  y ^)  profile: 

N( > - 1 )MJ  T - T ) 

e w el 

2(N+n(N*2)  (N*l)Te  | 
(3A) 

where  o*  is  found  from  the  incompressible  form  of 
the  momentum  equation  including  mass  transfer: 


Jo  _ (N*1)[N+2J  » 

7*  N I1 
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d-  * 
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and  x here  is  the  running  length  from  some  up- 
stream reference  point.  Thus  under  the  aforemen- 
tioned assumption  that  B is  a constant  over  some 
region  x)  - x ; X2  (and  zero  outside)  with  xj , x^ 
far  in  front  of  and  behind,  respectively,  the 
interaction  zone  x L,  and  using  Eq.  (2)  plus 
the  approximate  power  law  formula  Cf0  x -2 ( N+3 ) 
to  facilitate  analytical  integration,  Eq.  3B 
yields 


(N*3] 

1 • <V  B 

, 231 

:n+3! 

uKef/Te)Cjt 

''  21 

i"+2J 

n . 

(3C) 


which  in  conjunction  with  Eq.  3A  yields  J0. 

In  addition  to  tw,  mass  transfer  also  influ- 
ences the  profile  shape  away  from  the  wall  as 
given  by  the  following  turbulent  boundary  layer 
shear  stress  profile  relation  recommended  by 
Conrad  and  Donaldson4"  (similar  expressions  also 
have  been  proposed  by  others4^): 

7-  - 1 - 3n?  ♦ 3n3  ♦ (1  - nZ)(u/u  )(2B/C,  ) 

'w  0 

(O 

where  n y/S.  Further,  Eg.  4 is  to  be  used  with 
the  basic  turbulent  shear  stress  definition  that 
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Regarding  the  turbulent  kinematic  eddy  viscosity 
distribution  v?ff,  the  available  experimental 
evidence40-45  implies  that  its  functional  form  is 
significantly  affected  only  by  relatively  large 
surface  mass  transfer  rates  provided  the  mass 
transfer  effect  on  the  value  of  iw  Li  taken  into 
account;  to  a consistent  degree  of  approximation, 
then,  the  two-layer  piecewise-continuous  viscosity 
formulation  developed  for  use  in  interaction  prob- 
lems without  mass  transfer  by  lnger  and  Williams^' 
may  be  applied  also  to  the  weak-to-moderate  blow- 
ing or  surtion  cases  studied  here.  Thus  we  have 
for  i.  T : 
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where  r is  the  Kleb anoff  intermittency  factor  and 
f*  is  found  by  requiring  that  the  no  slip  con- 
dition Uo(0)  = 0 be  satisfied  upon  inward  inte- 
gration of  Eq.  5 from  the  outer  initial  condition 
Up ( 1 ) = Ue.  Thus  the  substitution  of  Eqs.  (4)  and 
(6)  plus  the  Crocco  integral  temperature  profile 
T(u)  into  Eq.  5 and  subsequent  integration  yields 
accurate  yet  fundamentally-based  incoming  turbu- 
lent boundary  layer  velocity  and  Mach  number  pro- 
files including  compressibility,  heat  transfer 
and  moderate  amounts  of  wall  suction  or  injection. 
The  results  satisfy  the  proper  boundary  conditions 
including  vanishing  gradients  at  the  boundary 
layer  edge,  conform  to  the  Law  of  the  Wall  near 
th?  surface,  are  continuous  across  the  entire 
boundary  layer  with  a velocity  defect-type  be- 
havior in  the  outer  part,  and  are  in  good  agreement 
with  experiment  over  a wide  range  of  transonic-to- 
moderately  supersonic  Mach  numbers^!  including  the 
effects  of  surface  mass  transfer^. 


3.  Numerical  Results  and  Discussion 
3. 1 Zero  Mass  Transfer 

To  provide  a basis  for  appreciating  and  scal- 
ing the  influence  of  mass  transfer  on  the  inter- 
action it  is  desireable  to  examine  first  some 
theoretical  results  for  the  impermeable  wall  case. 
Some  general  results  for  this  case  emphasizing 
the  Influence  of  Mach  number  have  already  been 


given  in  Ref.  11  and  thus  need  not  be  repeated;  we 
concentrate  here  on  more  recent  and  heretofore- 
unpublished  results  for  the  effect  of  Reynolds 
number  and  comparisons  with  experimental  data. 

The  predicted  influence  of  Reynolds  number 
on  the  interaction  pressure  field  for  a typical 
Mach  number  case  is  shown  in  Fig.  3.  It  is  seen 
that  there  is  a moderate  Reynolds  number  effect 
even  without  unseparation:  the  extent  of  the 
interaction  upstream  and  downstream  of  the  in- 
cident shock  decreases  with  increasing  Reynolds 
number,  tending  toward  a solution  typical  of  the 
response  to  a simple  step  pressure  rise  at  very 
high  Reynolds  number,  in  agreement  with  both  ex- 
perimental observations^.'?.  18,23-30  arid  Navier- 
Stokes  numerical  simulation^  of  turbulent  bound- 
ary layer-shock  wave  interactions.  Moreover,  at 
the  boundary  layer  edge  the  strengths  of  the 
local  shock  jump  and  post-shock  expansion  increase 
and  decrease,  respectively,  with  increasing 
Reynolds  number;  at  sufficiently  high  ReL  the 
post-shock  expansion  region  predicted  by  present 
theory  becomes  very  small  and  weak  and  hence 
probably  difficult  to  detect  experimentally. 

The  corresponding  upstream  influence  (defined 
as  the  distance  xup  ahead  of  the  shock  where  the 
local  interaction-induced  pressure  rise  is  only 
5%  of  the  overall  total)  at  various  shock  strengths 
as  a function  of  Reynolds  number  is  shown  in  Fig. 

4,  plotted  in  ratio  to  a0  (which  also  of  course 
experiences  Mach  and  Reynolds  effects).  These 
values  are  of  order  unity  (xup  — f0),  as  we 
should  indeed  expect  for  the  short-range  type 
of  interactions  characteristic  of  turbulent  bound- 
ary layers47(  and  decrease  markedly  with  both  the 
shock  strength  and  Reynolds  number,  at  moderate 
Reynolds  numbers.  xup/j0  decreases  monotonical ly 
with  Re^  approximately  as  a power  law.  It  is 
emphasized  that  in  addition  to  being  in  full 
qualitative  agreement  with  many  experimental  ob- 
servations, both  the  magnitude  and  parametric 
trends  of  these  theoretical  results  are  completely 
concordant  with  several  detailed  correlation 
studies'9-21 ,27  0f  upstream  influence  data  on 
interacting  turbulent  boundary  layers  that  directly 
verify  the  present  non-asymptotic  triple  deck 
flow  model . 

The  scale  effect  on  the  corresponding  inter- 
active displacement  thickness  growth  (Fig.  5)  is 
also  of  practical  interest  since  this  thickening 
often  has  a significant  back-effect  on  the  inviscid 
flow  and  shock  position  on  airfoils-*!  .48.50  or  in 
channel  flows** .37,49.  it  is  seen  that  the  pre- 
dicted displacement  growth  decreases  significantly 
with  increasing  Reyrolds  number  as  would  be  ex- 
pected (again,  this  trend  agrees  with  experiment). 

The  typical  distribution  of  interactive  skin 
friction  along  the  interaction  is  shown  in  Fig.  6 
and  illustrates  how  Cf  typically  decreases  to- 
ward the  shock  owing  to  the  adverse  pressure  gra- 
dient disturbance  induced  by  the  shock-boundary 
layer  interaction  (increasing  shock  Mach  number 
enhances  this  owing  to  the  stronger  local  inter- 
action pressure  gradient  involved").  When  the 
interaction  is  strong  enough,  the  present  theory 
predicts  vanishing  skin  friction  and  a very  short 
separation  bubble  slightly  behind  the  shock  foot, 
which  is  directly  confirmed  by  several  detailed 
studies  of  the  skin  friction  behavior  across  tran- 
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51.  Reynolds  number  has  the  expected  influence  on 
the  skin  friction  and  incipient  separation  behav- 
ior: the  relative  effect  of  the  interaction  at 
a given  shock  strength  decreases  with  ReL,  in- 
cipient separation  occuring  more  readily  at  lower 
Reynolds  number  as  observed  experimentally. 

To  further  validate  and  illustrate  the 
theory,  it  is  desireable  to  make  some  direct  com- 
parisons with  experiment  under  the  assumed  tran- 
sonic flow  conditions;  this  is  difficult,  how- 
ever, because  most  of  the  existing  transonic 
shock-boundary  layer  interaction  data  on  airfoils 
involve  high  local  Mach  numbers  (Mi  > 1.3-1. 4) 
with  a distinct  lambda-shock  interaction  pattern 
and  pronounced  boundary  layer  separation  which 
cannot  be  compared  meaningfully  with  the  present 
theory.  Nevertheless,  two  suitable  non-separated 
cases  were  found  in  some  published  NAE  wind  tun- 
nel tests  of  supercritical  airfoil  sections^2; 
the  measured  pressure  distributions  and  corres- 
ponding theoretical  predictions  (based  on  the 
local  pre-shock  Mach  number  and  Reynolds  number 
conditions  at  the  experimentally-observed  shock 
location)  are  shown  in  Fig.  7.  The  theory  is 
seen  to  predict  the  upstream  influence  well, 
whereas  it  overestimates  the  pressure  recovery 
downstream.  This  is  typical  of  such  airfoil 
tests  and  has  been  shown  31,48,50  to  be  caused  by 
the  fact  that,  in  contrast  to  the  normal  incident 
shock  theoretically  assumed,  the  actual  shock 
occurring  in  airfoil  experiments  is  usually 
oblique  (albeit  still  with  subsonic  post-shock 
flow)  owing  to  the  interactive  displacement  thick- 
ness back-effect  on  the  surrounding  inviscid 
flow;  this  lowers  the  actual  overall  shock  pres- 
sure rise  in  transonic  flow  20  - 30%  below  the 
normal  shock  value  at  the  same  incoming  flow  Mach 
number.  As  illustrated  by  the  good  comparison 
with  some  recent  DFVLR-Gottingen  interaction 
data50>51  on  a supercritical  wing  section  shown 
in  Fig.  8,  when  this  obliquity  effect  is  incor- 
porated the  present  theory  gives  a satisfactory 
account  of  the  interaction  downstream  as  well  as 
upstream  of  the  shock. 

Ackeret.  Feldman  and  Rott's  famed  experi- 
mental study! 2 of  shock-boundary  layer  inter- 
action on  a plate  and  wall  in  the  choked  transonic 
flow  of  a slightly-curved  wind  tunnel  nozzle  pro- 
vides some  further  examples  of  unseparated  tur- 
bulent flow  that  can  be  rendered  suitable  for 
comparison  with  the  present  theory;  we  have 
chosen  those  for  which  both  wall  and  inviscid 
pressure  distributions,  as  well  as  displacement 
thickness,  are  given.  It  should  be  emphasized, 
however,  that  direct  comparisons  with  their  data 
involve  numerous  uncertainties:  (a)  the  inviscid 
flow  edge  is  only  approximately  defined,  (b)  the 
shock  location  and  shape  are  uncertain  to  within 
.25  to  .50  x 60.  (c)  error  in  reading  the  curves, 
(d)  a significant  background  inviscid  pressure 
gradient  beclouds  interpretation  of  the  outer 
interaction  zone  and  the  incoming  turbulent  layer 
profile,  (e)  the  upstream  boundary  layer  history 
is  only  partially  understood,  especially  follow- 
ing forced  transition  cases,  and  (f)  a signifi- 
cant channel  blockage  effect  occurs  from  the  in- 
teractive boundary  layer  thickening,  which  re- 
duces the  effective  theoretical  shock  strength 
and  hence  the  downstream  interaction  pressure 
level3! ,37,49.  a typical  non-separating  inter- 


active pressure  field  measured  by  AF  & R12  is 
illustrated  in  Fig.  9a;  experimental  pressure 
distributions  along  both  the  surface  and  the 
approximate  boundary  layer  edge  are  compared  in 
Fig.  9b  with  our  theoretical  prediction  (cor- 
rected for  the  estimated  interactive  blockage  ef- 
fect using  the  Inger-Panaras  method^! >37) . while 
the  corresponding  displacement  thicknesses  are 
shown  in  Fig.  9c.  In  view  of  the  aforementioned 
uncertainties  from  the  data  and  limitations  of 
the  theory,  the  overall  agreement  is  seen  to  be 
quite  good.  In  particular,  the  following  defi- 
nitive features  of  the  interaction  theory  are  well- 
corroborated:  (1)  the  magnitude,  sign  and  stream- 
wise  extent  of  the  lateral  pressure  gradient 
effect  both  ahead  and  behind  the  shock;  (2)  the 
existence  of  a long  slow  interactive  pressure 
rise  (algebraic  rather  than  exponential ^ ) down- 
stream of  the  shock;  (3)  the  overall  streamwise 
scale  and  upstream  influence  distance;  (4)  the 
magnitude  and  shape  of  the  interactive  displace- 
ment thickness  growth;  (5)  the  local  inviscid 
pressure  jump  across  the  shock  at  the  boundary 
layer  edge;  (6)  a non-singular  inviscid  subsonic 
expansion  region  behind  the  shock  due  to  the 
viscous- inviscid  interaction. 

3.2  Mass  Transfer  Effects 
53 

Preliminary  study  indicated  that  the  domi- 
nant influence  of  mass  transfer  comes  from  the 
effect  on  the  profile  shape  away  from  the  wall; 
including  only  the  iw  - effect  gives  a significant 
error  in  both  magnitude  and  sign  of  the  profile 
changes.  The  typical  consequences  of  this  on  the 
interaction  solution  itself  are  illustrated  in 
Fig.  10,  which  shows  how  the  various  contributions 
to  the  suction/blowing  effect  on  the  Mach  no.  pro- 
file influence  the  wall  pressure  distribution 
(analogous  results  were  obtained  for  the  inter- 
active displacement  thickness  and  skin  friction). 
Whereas  the  contribution  of  the  mass  transfer 
effect  on  60  is  negligable  compared  to  that  on  ^w, 
the  influence  on  overall  profile  shape  is  large 
and  overwhelmingly  opposite  to  the  tw  - effect. 

This  conclusion  was  found  to  apply  over  a wide 
range  of  conditions  and  is  concordant  with  the 
fact  that  transonic  interactions  are  known  to  be 
very  sensitive  to  the  upstream  turbulent  boundary 
layer  profile  form  factor37. 

Referring  hereafter  to  the  complete  mass 
transfer  model  we  observe  that  suction,  because 
of  its  predominant  effect  in  decreasing  the  Mach 
number  gradient  and  hence  enhancing  the  profile 
"fullness"  away  from  the  wall,  reduces  the 
streamwise  extent  and  thickening  of  the  inter- 
action, making  it  appear  more  inviscid-like  in 
character  with  a steeper  adverse  pressure  gra- 
dient. Thus  suction  is  qualitatively  equivalent 
to  an  increase  in  Reynolds  number.  Blowing  has 
the  expected  opposite  effects  of  spreading  out 
tne  interaction  pressure  field  and  increasing  the 
displacement  thickness. 

Results  of  a systematic  study  of  suction/ 
blowing  effects  on  the  interaction  pressure  field 
for  a typical  case  are  presented  in  Fig.  11.  In 
Fig.  1 1 A it  is  seen  that  moderate  amounts  of  suc- 
tion significantly  reduce  the  upstream  influence 
distance  and  overall  streamwise  extent  of  the 
interaction  and  steepen  the  adverse  wall  pressure 
gradient,  whereas  blowing  has  equally  the  opposite 
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effect.  Concordant  with  these  trends,  suction 
also  strengthens  *he local  shock  jump  at  the 
boundary  layer  edge  while  reducing  (perhaps  even 
eliminating  at  high  enough  B)  the  degree  and  ex- 
tent of  the  post-shock  expansion  region  (Fig.  1 1 E ). 
The  non-dimensional  upstream  influence  distance 
versus  the  mass  transfer  parameter  B is  plotted 
in  Fig.  12. 

The  predicted  influence  of  suction  and  blow- 
ing on  the  interactive  displacement  thickness 
distribution  is  shown  in  Fig.  13,  while  the  total 
downstream  thickness  is  plotted  vs.  B in  Fig.  14. 
It  is  seen  (as  expected)  that  this  thickening  is 
significantly  influenced  by  mass  transfer;  for 
example,  the  moderate  suction  value  B = -0.0003 
reduces  Ai*(«)/s0*  nearly  three-fold.  It  is 
noted  that  the  influence  of  Reynolds  number  on 
the  interaction  in  the  presence  of  surface  mass 
transfer  effects  described  below  was  also 
studied^;  suffice  it  here  to  state  that  over 
the  range  of  values  |B|  < .0005,  its  relative 
effect  was  found  to  be  quite  similar  to  that 
shown  above  in  Section  3.1. 

The  interactive  local  skin  friction  is  of 
particular  interest  since  it  identifies  possible 
flow  separation  and  provides  re-initialization 
data  for  continuing  boundary  layer  calculations 
downstream  of  the  interaction  zone.  We  note  in 
this  regard  that  although  the  present  theory  is 
no  longer  valid  for  separated  flow  [Cf ( x ) <•  0 
over  some  portion  of  the  wall]  it  is  still  useful 
to  indicate  trends  toward  this  situation,  i.e., 
where  and  when  incipient  separation  (Cf  - 0 at 
some  x)  first  occurs.  As  indicated  in  Fig.  15, 
the  influence  of  mass  transfer  involves  two  op- 
posing effects;  far  upstream  or  downstream  the 
skin  friction-increasing  effect  of  suction  domin- 
ates (which  tends  to  delay  separation)^ , whereas 
in  the  neighborhood  of  the  shock  foot  (|x/a0|<  1) 
the  suction-induced  steepening  of  the  local  adverse 
pressure  gradient  becomes  of  controlling  impor- 
tance and  Cf  is  actually  reduced.  Thus  in  con- 
trast to  what  occurs  in  non-interacting  boundary 
layers,  slight  suction  here  actually  has  a locally- 
adverse  effect  in  hastening  interactive  incipient 
separation  under  the  shock;  rather,  it  is  small 
amounts  of  blowing  which  can  delay  (and  in  the 
present  example  eliminate  altogether)  such  separa- 
tion. When  the  blowing  is  sufficiently  strong. 
Fig.  15  indicates  that  the  upstream  skin  friction 
reduction  will  eventually  exceed  that  in  the 
shock  foot  region,  whereupon  Cf  • and  hence  the 
incipient  separation  point  begin  to  move  forward 
in  front  of  the  shock.  Presumably,  at  even  high- 
er blowing  rates  separation  ultimately  will  occur 
well  upstream  of  the  shock;  however,  the  present 
theory  must  be  modified  to  study  this  question 
since  several  of  the  underlying  simplifying  as- 
sumptions become  doubtful  at  these  rates  (BRIO'S). 


layer  normal  shock  interactions1 ’^).  Thepresent 
theory  indeed  shows  a small  separation  region  ap- 
proximately one  boundary  layer  thickness  in  length 
under  the  shock  foot,  in  rough  agreement  with  ex- 
perimental findings.  Moreover,  it  is  seen  that 
this  can  be  eliminated  by  small  blowing  rates  ex- 
ceeding B • 5 x 10*3  fo  produce  a fully  attached 
unseparated  flow  throughout  the  interaction  zone, 
albeit  one  which  emerges  downstream  with  an  ex- 
tremely low  Cf  level  that  is  completely  out  of 
local  equilibrium. 


It  is  noted  that  while  there  presently  exists 
no  systematic  quantitative  experimental  data  on 
transonic  shock-boundary  layer  interactions  with 
surface  mass  transfer,  some  observations  in  a 
channel  flow  test  have  been  reported^  that 
qualitatively  support  the  foregoing  conclusions. 
Some 

compressi 

also  observe  our  theoretically-predicted  reduction 
of  interactive  thickening  and  upstream  influence. 
However,  to  the  authors'  knowledge,  there  is  no 
data  regarding  the  influence  of  mass  transfer  in 
the  shock  foot  region  of  transonic  shock-unsepara- 
ted  turbulent  boundary  layer  interactions. 


recent  experimental  results  on  a supersonic 
ressive  interaction  flow  with  wall  suction” 


4.  Concluding  Remarks 

This  non-asymptotic  study  of  weak  normal 
shock-turbulent  layer  interactions  for  two-dimen- 
sional non-separating  flows  including  mass  trans- 
fer has  shown  that  even  small  amounts  of  suction 
(-mw/Peue  ; 5 x 10*4)  appreciably  reduce  both  the 
streamwise  scale  and  thickening  effect  of  the 
interaction  but  hasten  the  onset  of  separation 
slightly  behind  the  shock  foot.  Equal  amounts  of 
weak  blowing  on  the  other  hand  can  completely 
eliminate  interaction-induced  separation.  These 
results  were  found  to  be  mainly  a consequence  of 
the  mass  transfer  effect  on  the  incoming  boundary 
layer  Mach  number  profile  shape  which  in  turn 
significantly  affects  the  interaction  pressure 
distribution  and  hence  the  local  skin  friction. 

The  present  theory  provides  a useful  analytical 
framework  for  the  evaluation  and  parametric  study 
of  the  interaction  mass  transfer  effects  in  a 
variety  of  practical  applications.  Moreover,  it 
provides  the  basis  for  further  improvement:  ex- 
tension of  the  analytical  model  to  larger  mass 
transfer  rates  by  including  their  explicit  effect 
on  the  viscous  disturbance  sublayer  and  mean  tur- 
bulent boundary  layer  eddy  viscosity  equations. 
With  such  added  features,  it  should  be  possible  to 
examine  the  basic  question  of  how  mass  transfer 
influences  incipient  separation  over  a very  wide 
range  of  suction  or  blowing  rates  as  well  as  Mach 
and  Reynolds  numbers  without  the  need  for  present- 
day  empiricsms. 


The  aforementioned  favorable  local  effect  of  Acknowledgment 
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Thus  in  the  large  scale,  suction  always  has  the 
beneficial  influence  of  promoting  a more  rapid 
equilibration  of  the  turbulent  skin  friction 
downstream. 


1 It  should  be  noted  that  the  often- important 
blockage  corrections  in  such  tests  due  to  inter- 
active thickening  of  the  wall  boundary  layers  can 
be  significantly  reduced  by  suction  (see  Fig.  13). 
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